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Abstract 

We formulate the theory of the two-stream instability (e-cloud insta- 
bility) with electrons trapped in quadrupole magnets. We show that a 
linear instability theory can be sensibly formulated and analyzed. The 
growth rates are considerably smaller than the linear growth rates for 
the two-stream instability in drift spaces and are close to those actually 
observed. 

1 Introduction 

The Proton Storage Ring (PSR) at Los Alamos has been troubled for some time, 
P] , [2], [2], 0], [5], [6] , [7], [8] , [9], by an instability that is probably a two-stream 
instability of the proton beam with background electrons, i.e. an electron-cloud 
instability. We have previously considered the possibility that the instability 
for a bunched beam occurred because of electrons in drift regions that were 
renewed from turn to turn, [10) ; in this case the phase memory of the coherent 
motion has to reside in the proton beam which excited the fresh electrons on 
each turn which then drove the tail of the proton bunch to larger amplitudes. 
In this note we will consider instead the possibility that the instability is due to 
electrons that survive from turn to turn. The most likely place in the ring where 
the electrons can survive with coherent phase information from turn to turn is 
in the quadrupoles, where they are trapped in the magnetic mirrors formed 
by the cusp-shaped fields and can drive the e-p instability in a similar way to 
free electrons. In this note we will present a simple model of this two-stream 
instability with electrons trapped in quadrupoles. 



"This work was supported by the US Department of Energy under Contract Number DE- 
AC52-06NA25396. 
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1.1 Electron Trapping and Dynamics 



A major assumption of this note is that there are abundant electrons in the PSR; 
experimentally this has been observed, though the source is not completely 
clear. It is likely that some form of beam induced multipactor gives rise to 
the electrons, perhaps initiated by a very small number of lost beam particles, 
though other explanations are possible. Normally, one would expect that with 
a bunched beam electrons would be expelled during the beam gap and that 
one could not have an e-p instability; however, the electrons, however they are 
produced, cannot be driven quickly to the walls in the quadrupoles which act 
in the transverse direction as very effective magnetic mirrors. It thus seems 
possible that electrons in the quadrupoles could drive the e-p instability. To 
investigate this possibility further in this section we will make simple estimates 
of the electron motion in quadrupoles to establish that electrons can be trapped 
there for multiple turns and thus carry coherent phase information to drive the 
instabillity. A more accurate investigation of the electron motion in the complex 
geometry can and should be done using computer codes, [TT] . 

The dominant aspect of electron motion in the quads is the rapid rotation 
about the magnetic field lines; the cyclotron frequency is 

'-555? (1) 

where e is the charge, B is the magnetic field, m is the mass, and c is the speed 
of light. For electrons we have 



f c = 2.8B GHz, (2) 

where B is the magnetic field measured in kilogauss. Thus, even very low fields 
near the axis give rise to cyclotron frequencies that are hundreds of MHz; most 
electrons will have cyclotron frequencies that are multiple GHz. The radius of 
this rotational motion is, for electrons, 

p = 3.37* !0~ 3 ^- cm, (3) 

where £ is the transverse electron energy in eV. Only very energetic electrons in 
low field regions will have gyroradii approaching 1 cm; most will have gyroradii 
that are much less than 1 mm. Electrons are thus confined transversely to 
the magnetic field on cyclotron orbits of small radii and many are confined 
longitudinally (for electrons) along the magnetic field by the increasing magnetic 
field with radius, i.e. by 'mirror' confinement. (Note that longitudinal for the 
electrons is transverse to the beam direction.) 

Of course, particles with large components of velocity parallel to the mag- 
netic field, i.e. those in the 'loss cone', are not confined; presumably these give 
rise to the electron 'tracking' that has been observed in the quadrupoles. We 
will ignore the rapid electron cyclotron motion in the quads and concentrate 
on the longitudinal electron mirror motion and transverse drifts due to electric 
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fields and to magnetic field non-uniformity, i.e. a 'guiding center' description of 
the trapped electrons. 

In the transverse direction (for electrons) there are three components of elec- 
tron drift, that due to the gradient in the magnetic field, the so-called VB drift, 
that due to the field line curvature, and that due to any electric fields that are 
present, the E x B drift. These drifts give rise to electron velocities perpendic- 
ular to the magnetic field; in fact, in the quads, the drifts are along the direction 
of the beam axis and thus can lead to electron loss out the ends of the quads. 

The V-B drift and curvature velocities are given by 

+ 2v f\) b x VB 
= 2eB B C ' (4) 

where & is a unit vector in the direction of the magnetic field, the v\ term is 
due to the gradient drift and the term is due to the curvature. If we assume 
the parallel and perpendicular electron velocities to be roughly the same and 
adopt the usual model of quadrupolc magnetic fields in which a component is 
linear in transverse displacement from the axis, i.e. 

B = B'r, (5) 

then, defining the WB confinement time, Tvs to be the time for an electron to 
drift half the length of a quad, Lq , we get 

_ eB'r^LQ 

TVB 4~Ec ' (6) 
where E is the thermal energy of the electron. This becomes 

Tvb = Br ^ Q Msec, (7) 

where B' is the field gradient in T/m, f is the radius in cm, £ is the energy in 
eV, and Lq is the quad length in cm. As an example typical of the PSR, if we 
take B' = 3.7, f = 2.5, and Lq = 47, then 

272 

Twb = — Msec. (8) 

Note that this is an overestimate of the drifts since the actual drift reverses sign 
as the electrons move out along the magnetic field lines toward the poles. If the 
electrons only have energies that are a few hundred eV then the confinement 
time is tens to hundreds of turns and is probably longer than the growth time 
for the e-p instability. 

The E x B drift velocity is given by 

V EXB = (9) 

The electric field is due to the proton beam and to any electrons that are present. 
The electric potential due to the proton beam alone is given by 
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where /? is the beam velocity scaled by the speed of light and / is the (time- 
dependent) beam current. The beam current varies by 100% in one revolution 
period (the beam is bunched), but we will estimate drifts using the average 
current and resulting field. Note that electrons spend a lot of time near the 
magnetic mirror points where we expect that the E field will mostly be parallel 
to the B field and will give rise to only small drifts. Nevertheless, the E x B drift 
velocity due to this term alone, assuming it acts all the time, would give an 
electron confinement time of 

T ExB = 8.34 L Q nsec, (11) 

where the current, T, is measured in amps, if we again take B — 3.7, f = 2.5, 
T = 10, and Lq = 47, then TexB = 761 nsec, i.e. electrons would be confined 
for several turns, even with this overestimate of the E x B drift. With a more 
realistic calculation, including the full orbit dynamics of the electrons and the 
reverse drifts that occur when only the electrons are present, it is likely that the 
electrons will be confined for many turns. 

Electrons to the left and right of the beam, horizontally, are free to move 
vertically (initially) until they move out radially along the field line to a region 
of greater field strength. Electrons above and below the beam, vertically, are 
free to move horizontally (initially) until they move out radially along a field line 
to a region of greater field strength. A complete model of the electron motion 
is very complicated, but a simple model will suffice to treat the motion of the 
center of mass of the electrons for oscillations near the beam axis. Let us note 
that for electrons that can move vertically, i.e. those to the left and right of the 
beam, the restoring mirror force exactly vanishes at zero vertical position and 
the restoring force reverses sign there. For electrons that can move horizontally, 
i.e. those above and below the beam, the restoring mirror force exactly vanishes 
at zero horizontal position and the restoring force reverses sign there. Thus, 
in both transverse directions we should expect the restoring potential for an 
electron to be approximately a harmonic oscillator potential near the axis. To 
see this in more detail, let us begin with the equation from Krall and Trivclpiece, 
[12] , for the equation of motion along a field line of a particle in a magnetic field 

dt 2 ~ 2B ds ' 1 ' 

where s is the distance along the field line, v±q is the initial value of the trans- 
verse velocity, and Bq is the initial value of the magnitude of the magnetic field. 
The components of the quadrupole field are 



B x = B' y, 



(13) 



B y = B' x. (14) 

We thus see that 

dB = 2B' xy _ 
ds x 2 + y 2 

From this we see that a particle that starts at x = xq, y = satisfies the 
approximate equation 



fy M { b'qv 2 

dt 2 BqXq 



^«-(^)». ( 16 ) 



i.e. it is approximately a harmonic oscillator with a squared angular frequency 
of 



' „,2 



2 _ ^O^iO 



— 



BqXq 

But -Bq ~ B' xq, so 



(17) 
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It thus appears that modeling the mirror trapping of the electrons by a harmonic 
oscillator potential, but with a large spread in oscillation frequencies, should be 
a fairly good approximation. 



2 Dipole Model of the e-p Instability 

In this section, in order to find thresholds and growth rates, we will present 
a simple theory of the e-p instability. The model for the linear theory of the 
instability in this section that we use is similar to the theory of Keil and Zottcr, 
[13] . We model the proton beam by the beam centroid at each azimuthal po- 
sition around the ring. The background electrons have a complex distribution 
both in physical and in velocity space determined by their formation, capture 
in the quadrupoles, interaction with the proton beam, and loss, as discussed in 
the previous section. We cannot hope to accurately model all of these effects 
in an analytic theory; we will simply assume that the electrons have a distri- 
bution in the squared magnetic bounce frequency, g m — uj 2 ^, and that at each 
bounce frequency those electrons are described by their centroid position, with 
electrons at a different bounce frequency having a different centroid. We assume 
the proton beam moves at a constant azimuthal velocity around the ring and 
is subject to a constant transverse focusing force that produces betatron oscil- 
lations at the betatron frequency, i.e. we make the smooth approximation, [14] , 
We only model proton beam and electron motion in one transverse direction. 
The protons and electrons are assumed to interact with each other via a force 
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that is linear in the relative displacement of the centroids of the protons and 
electrons. The equations of motion for the centroids are thus given by 

d d d d 



dt 2 



+ u m y em = uj e (y p - y em ) (20) 



where y p {9,t) is the proton ccntroid position at an azimuth, 9, around the 
machine and time, t, lu is the proton beam angular revolution frequency in the 
machine, and u>p is the angular betatron frequency of the protons. The proton 
beam centroid only responds to the net electron centroid position, Y e , which is 
given by 



f(9m)yem(0,t)dg m , (21) 



where f(g m ) is the equilibrium distribution function of electrons in the squared 
bounce frequency and y em (6,t) is the centroid of electrons with a particular 
bounce frequency. The coupling frequencies ui p and u> e are given by 

, = 2iy e c 2 

e nb(a + b)R K ' 

4 . (^H (23) 

with N p the number of protons in the machine, r e the classical electron radius, c 
the velocity of light, 7 the relativistic factor of the proton beam, a and b the sizes 
of the proton beam, F the neutralization fraction of electrons, and R the effective 
radius of the ring. Note that the inter-species force is assumed to depend linearly 
on the distance between the beam centroids; this is approximately correct for 
small amplitudes of oscillation, but clearly fails at larger oscillation amplitudes. 

Also note that we have inserted a linear damping term with coefficient Yd into 
the proton equation to account for the chromatic spread in proton revolution 
frequencies; the different revolution frequencies will give different longitudinal 
velocities which will Landau damp the transverse oscillations. A more extensive 
model would have the proton beam described by a distribution function in 
the azimuthal direction and take into account the Landau damping due to the 
spread in azimuthal velocities. The approximation we have adopted mimics 
this damping and has the same functional dependence as the result of this more 
extensive model (see below), i.e. the damping depends on 1) the energy spread, 
2) the momentum compaction factor, and 3) the mode number (through the 
derivative in the damping term). Thus, this damping term will give rise to 
the correct qualitative behavior with the correct functional dependencies, i.e. 
damping of off-axis oscillations as they phase-mix away. We can estimate this 
damping rate of transverse oscillations due to this spread to be the chromatic 
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fractional tune spread times the betatron frequency. Note that the chromatic 
fractional tune spread is just the chromaticity times the energy spread, i.e. it 
measures the longitudinal velocity spread and its influence on the transverse 
oscillations. We do not include the transverse tune spread due to space charge 
and machine nonlinearities because we are using a dipole model and the centroid 
motion of the protons does not depend on these terms. 

IW£)£ (24, 

Because we are using an unbunched beam model, i.e. the smooth approximation, 
the average neutralization around the ring will be smaller than the neutralization 
in the quadrupoles by roughly the ratio of the ratio of total quadrupole length 
to the ring circumference; thus the neutralization fraction in a quadrupole will 
be about 20 times F since quadrupoles are about 10% of the circumference and 
only about half the electrons can move vertically. 

We have seen in the context of the drift space instability model, ,10], that 
bunching doesn't have a large effect on the instability, and we assume the same 
to be true here. There seems to be no simple way to incorporate bunching; a 
moderately realistic model would result in a dispersion equation which would 
be an infinite matrix eqation with all unbunched beam modes coupled. The 
unbunched beam model of this paper would then be just the diagonal approxi- 
mation to this matrix equation. It is likely that an extensive numerical investi- 
gation would be required to resolve the behavior. 

The above model is overly simplified, but contains most of the important 
physics. It will break down, of course, if the electron loss rate is too high. Of 
couse, we are also assuming that the background electron density, on average, 
is constant so that if electron generation and loss rates fluctuate rapidly our 
model should fail. 

The model of Bosch, [15., for the effect of beam gaps on the trapped ion 
instability in an electron ring also considers the effect of a large spread of (ion) 
oscillation frequencies on the instability, and his formulation is similar to ours. 

If we assume that the perturbations have a dependence on time and angle 
proportional to e - lut + m0 j then the equations become 

(-(u - nuj ) 2 - iT d (uj - nuj Q ) + u 2 + u 2 v )y v = u 2 v Y e (25) 

(^e + W m - w2 )^ = ^lv P - (26) 

Solving equation [26] for y e and using equations [21] and [25] we find 



((uj - nw ) 2 + iT d (u - nuj a ) - (Jp - lu 2 ) = -uj 2 uj 2 f ^} 9 ™^ j^ffm, (27) 

where we have used the definition of g m — u^. We have to deal with the 
singularity in the integral on the right hand side of this equation. We adopt 
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the Landau prescription, see [12) , where the integral is replaced by the principal 
value plus iri times the residue at the pole; 



/ f (g 2 n) 2 d 9m = / ^ 2 + ^ ~ (28) 

J g m + ujj-uj 2 J g m + ujj-uj I 

Equation [23 thus becomes 

((lu - nuj Q ) 2 + iT d {uJ - ntJo) — Up — lu 2 ) 



where the bar through the integral sign indicates principal value. This is the 
dispersion relation for the two-stream mode. To solve it we have to specify 
the distribution function of electron bounce frequencies, /. Of course, there 
should be no electrons in the 'loss-cone', i.e. at zero w^, but otherwise the 
detailed distribution depends on their formation, capture in the quadrupoles, 
interaction with the proton beam, and loss. We will simply take one distribution 
as an example, one in which the distribution is constant between a minimum 
squared bounce frequency and a maximum squared bounce frequency; i.e. 

/(<?m) = 9min — Qrii — Smax-j 

9max 9min 

otherwise (30) 

With this distribution the dispersion equation, [2TJ[ becomes 

(u — nuj ) 2 — u 2 +u 2 — ir d (u) — nuj ) 

Ue^p ,gmax-U 2 +UJ 2 
— (H 2 _ 2 _ )+7Tl) (31) 

Though this is a transcendental equation and can't be solved exactly, we note 
that the coefficient of the logarithmic term is small and the logarithm varies 
slowly, so we can simply solve iteratively. The remainder of the equation is a 
quadratic for u> — uluq and the resulting approximate solution is 

uj ncj — 

^^(ln( f--^ 2 o)2+ ^ 2 ) + * (32) 

Qmax Qmin {^^0 ) LU e Qvain 

We note the damping due to the phase mixing term, as expected, and the 
usual upper and lower sidebands. Note that we have taken oj ss udq inside the 



2 2 f Agm) , 
-mululf{u 2 -u% (29) 
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logarithm because the mode numbers are usually rather high 30 — 50 and this 
is a good (few percent) approximation for the real part of the frequency. Let us 
expand just the imaginary term under the square root to find the damping and 
growth rates. For convenience define the real frequency shift to be 



LJ 



1 \aI 2 , 2n -n2 4w e w p , 9max ~ (rUJo) 



[A(^ + oj p )-T d -- — — (ln( ))] (33) 

^ ymax ymm \ lt ^0) w e yra%n 



Note that to a good approximation lo s w w^. Expanding the imaginary term in 
the square root we get 



U) re ncj - -z- ± w s (l - — jt- — rr) (34) 

\9max ymin ) ) 



Note that the upper side band (plus sign) is always damped, but that the lower 
side band can be unstable if 

2 2 

777 ^7 ^ > T d , (35) 

\9rnax Qrain ) ) 

with growth rate given by 

Tgrowth = TT—, 7 7: - "T"- (36) 

A number of modes in lower side bands can be unstable, limited only by the 
condition oj 2 + g m in < uj 2 < ojI + g ma xi with roughly equal growth rates (there 
is some weak dependence on mode number in u> s ) and this is consistent with 
experiments where multiple modes are usually observered, |16| . 

2.1 Example 

Let us look at an example typical of the PSR; let us take 

a = b = 1.8 cm, 

lu = 2tt * 2.8 MHz, 
ujp = 2.2 * u/q, 
2irR = 89 m 

If we express the number of particles in the ring as 



then we can compute 



iV p = Wx 10 13 , (37) 



2 = 1.76* A/"* 10 17 sec" 2 , (38) 
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and 



lj 2 p = 0.5182 *F*Af* 10 14 sec" 2 . (39) 

In the PSR the measured vertical chromaticity is about —1.68 and the energy 
spread (typical conditions) is about 0.5% so we take the chromatic tune spread 
to be about 0.009, i.e. a fractional tune spread of 0.43%, then 

r d w 0.0252* 10 6 sec" 1 . (40) 
We take the frequency shift to be 

w s «cj /3 = 3.9*10 7 sec" 1 . (41) 

To estimate g m ax we use equation [18J setting the maximum transverse energy 
to the beam potential; the result is 

g max « 3.48 * N * 10 17 sec" 2 , (42) 

where we used Xo ~ a = 1.8 cm. Note that we simply ignore g m in, i-e. assume 
it is zero; it only modifies our results by a small factor. 

If we evaluate the threshold condition, equation [331 using equations |3"51 13TJ1 
BUI HU and 221 we find the criterion for instability to be 

F*N > 0.188; (43) 

in other words, once the product of the particle number (times 10 13 ) and per- 
cent neutralization is about 19.0, we can expect instability. Recall that the 
neutralization fraction in quadrupoles will be about 20 times higher than F 
since quadrupoles are only about 10% of the ring and only half the electrons 
can move vertically. At threshold the growth time is infinite, but if, for simplic- 
ity, we assume that we are a factor of 2 above the threshold, F*Af — 9.4* 10~ 2 , 
then using [Ml [Ml HOI SD and gj in equation [Ml we find 

7growth « 12.6 KHz, (44) 

i.e. a growth time of about 222 turns. These estimates are only intended to show 
that the results seem to be within a factor of two or three of the observations 
and that the theory is thus a possible explanation of the observed instability. 



3 Discussion 

Our results show that electrons trapped in quadrupoles are a plausible explana- 
tion of the two-stream instability observed in the PSR. The growth times found 
are considerably closer to the observed values than the linear growth times de- 
rived from the instability treatment for electrons in drift spaces, QJj]. The reason 
for this is that the electrons confined in quadrupoles have a very large frequency 
spread due to the wide variation in magnetic bounce frequencies as compared 
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to those in drift spaces which have only a very small spread in space charge 
confinement bounce frequencies. Thus, many fewer electrons are resonant at a 
particular frequency. 

In addition, if the instability is due to electrons trapped in quadrupoles, 
then the transverse momentum kick given to the protons is easily explained; the 
momentum is transferred from the quadrupoles via the electrons, rather than 
having to be transferred only from electrons, as in the drift space theory. 

Clearly a great deal more work can be done to refine this model. A kinetic 
description of the proton beam could be used, and would give a more sensitive 
dependence of the phase-mixing damping that depends on the detailed pro- 
ton distribution. An investigation of different electron distribution functions, 
perhaps motivated by detailed simulation of electron formation and capture dy- 
namics, would give threshold and growth rate estimates that are better founded 
than those in this note. The formulation of a bunched beam model would 
be considerably more difficult, but might be worthwhile. Finally, a composite 
model with both drift space electrons and quadrupole trapped electrons would 
be very difficult to analyze but might be necessary to fit all observations in real 
machines. 
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